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ABSTRACT 

Context. It is generally agreed that hydrogenation reactions dominate chemistry on grain surfaces in cold, dense molecular cores, 
saturating the molecules present in ice mantles. 

Aims. We present a study of the low temperature reactivity of solid phase isocyanic acid (HNCO) with hydrogen atoms, with the aim 
of elucidating its reaction network. 

Methods. Fourier transform infrared spectroscopy and mass spectrometry were employed to follow the evolution of pure HNCO ice 
during bombardment with H atoms. Both multilayer and monolayer regimes were investigated. 

Results. The hydrogenation of HNCO does not produce detectable amounts of formamide (NH 2 CHO) as the major product. 
Experiments using deuterium reveal that deuteration of solid HNCO occurs rapidly, probably via cyclic reaction paths regenerat¬ 
ing HNCO. Chemical desorption during these reaction cycles leads to loss of HNCO from the surface. 

Conclusions. It is unlikely that significant quantities of NH 2 CHO form from HNCO. In dense regions, however, deuteration of HNCO 
will occur. HNCO and DNCO will be introduced into the gas phase, even at low temperatures, as a result of chemical desorption. 
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1. Introduction 

Isocyanic acid, HNCO, is the simplest molecule containing 
the four most abundant elements: hydrogen, carbon, nitro¬ 
gen, and oxygen. It is an important interstellar molecule, 
with gas phase abundances of around 10 “® - 10 “^ with re¬ 
spect to H 2 in m olecular clouds, wh ere it is believed to trace 
dense, cold gas (iJackson et al.lll9^ . S ince its first detectio n 
in the Sgr B molecular cloud complex (ISnvder & Buhll[r972h . 
HNCO has been detected in multiple environme nts, includ¬ 
ing hot cores (e.g. iHelm ich & v an Dishoecldll997l) . high mass 
young stellar objects (iBisschop et al.l2007a ) , molecular out flows 

(IRodrfguez-Fer nandez et al.||.2O10|), comets dUis et al nm, and 

other galaxies dNguven-O-Rieu et al1ll991l) . It has been shown 
that HNCO is also a tracer of warm gas, and that its forma- 
tion origins are likely to be predominantly in icy grain mantles 
(iBisschop et al.ll^07al) . 

Isocyanic acid has not yet been observed in the solid state, 
but it is believed to be responsi ble for the formation of the abun¬ 
dant c yanate ion OCN~ (e.g. ISoi fer et al.l 19791: iDemvk et ^ 
119981: iLowenthal et"^ 120021 1 van Broekhuizen et al 
via reaction of the acidic HNCO with bases 


120051) 
such 


as 


NH^ (Raunier et al. l2003al Ivan Broekhuizen et a 

.1 12004 

Mispelaer et al.l 12012l) 

and HiO (Raunier et al. 

l2003b 

Theule et al. 1201 lb, or ffi 

/ the irradiation of ices with ultra- 

violet photons (jEacy et al.||l984b or protons (iMoore et al.lll983h. 


Moreover, its solid phase chemical reactions have been shown 
to give rise to isomerisation, with cyanic acid, HOCN, form ed 
thermally in mixtures of HNCO and H 2 O dTheule et al.ll201 ll) . 


The chemical network surrounding HNCO has not been fully 
studied experimentally. Early theoretical m odels assumed that 
HNCO formed only in the gas phase (e.g. Ilglesiaslll977h . but 
more recent studies contend that HNCO forms on grain surfaces 
via the the rmal reaction NH + CO or by the hydrog enation of 
OCN (e.g. iGarrod et al.ir2008t iTideswell et ^fcoKlI) . Its pres¬ 
ence in the gas phase is explained by the subsequent desorption 
of HNCO from grains, or by the destruction of larger species 
such as urea, (NH 2 ) 2 CO. Abundances of HNCO are enhanced 
in shocked regions, and thus it is likely that sputtering or addi- 
tional gas phase forma tion routes are active in post-shock gases 
(IZinchenko et al.ir 200 (]l) . 

The abundance of solid HNCO is predicted to remain rel¬ 
atively low, at ~ 10 with respect to H 2 O, and thus is not 
detectable in infrared spectra. If the d ominant formation route 
of HNCO (NH -H CO) were efficient (iGarrod et al.ll2008l) . de¬ 
struction routes would be required to explain this. As mentioned 
above, the formation of OCN“ by reaction of HNCO with H 2 O 
and NH 3 are favourable routes, the reacti ons having activatio n 
ban-iers of 26 + 2 kJ moU' (3130 + 240 K:lTheule et al.l(l201 ll) ) 
and 0.4 ± 0.1 kJmoF* (48 + 12 K: iMispelaer et al.1 (l2012l) ).'re- 
spectively. Irradiation of HNCO with vacuum ultraviolet (VUV) 
radiation has been shown to produce formalde hyde (H 2 CO), 
formamide (NH 2 CHO), and urea (H 2 NCONH 2 ) (iRaunier et al.l 
12001 . 

In dense molecular clouds the secondary photon field is weak 
and hydrogen atoms have a long resid ence time on grain sur¬ 
faces because of th e low temperature (ITielens & HagenI [19821: 
lAmiaud et al.ll200l . Hydrogen atoms are mobile on the surface 
at 10 K, and thus hydrogenation reactions dominate low temper- 
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ature ice chemistry. Experimentally, there still remain many hy¬ 
drogenation reactions to characterise, although studies have been 
carried out, particularly on simple molecules. The hydrogenation 
of C O to form H 7 CO and CH^OH has been extensively studied 
(e.g. Hiraoka et alJll994l : IWatanabe & Kouchill2002l:iFuchs et all 


[ 20 ^ as have the re action pathways to the formation of H 7 O 


from atomic oxyge n dHiraoka et al 


molecular oxyge n ( Mivauchi et al.l 200^ IoppoIo et all l2008t 
IChaabouni et alj|2012l) . and ozone ( Mokrane et al. 200^ . The 


|199M Dulieuet^ 


formation of C7H5 OH, CH4, H^CO and CH c,OH from CH3CHO 
has been studied bv iBisschop et alJ (l2007bll. while the forma tion 
of NH2OH by the hydrogenation of NO dCongiu et al.l2012ll and 
the fo rmation of CH2NH and CH3NH2 from HCN ( Theule et alJ 
I 2 OI Ih have recently been demonstrated. 

Given that atomic hydrogen is present at fractional abun¬ 
dance s of [Hi]/[H 2 ] ~ 10“^ in molecular clouds (iLi & GoldsmithI 
l2003h . hydrogenation is likely to dominate the destruction path¬ 
ways of HNCO at low temperatures. Theoretical studies may 
suggest that the radical intermediates HNCHO or NH 2 CO 
form rapidly from HNCO, and the stable molecule formamide 
forms on further h ydrogenation of these intermediate species 
dGarrod et al.ll20d^ . giving overall 


HNCO + 2H ^ H 2 NCHO. (1) 

This work focuses on the reaction of HNCO with H, which 
has not yet been examined experimentally. The experiments per¬ 
formed are introduced in §| 2 ] the results of these experiments are 
presented and discussed in § [3] while the astrophysical implica¬ 
tions of the results are considered in §| 4 l 


2. Experimental 


Experiments were performed using two differen t experimen¬ 
tal set-ups: RING, as describ ed in TheuleeLak (1201 ih . and 
EORMOLISM, as described in lAmiaud et al.l (l2006l) . 

The RING set-up was used to perform multilayer, bulk 
ice experiments. Briefly, RING consists of a gold-plated cop¬ 
per surface within a high vacuum chamber (a few 10 “® mbar). 
Molecular species in the form of room temperature gas are dosed 
onto the gold surface (15 - 300 K) by spraying via an injection 
line. The infrared spectra of the molecular solids are recorded by 
means of Eourier transform reflection absorption infrared spec¬ 
troscopy (FT-RAIRS) using a MCT detector in a Vertex 70 spec¬ 
trometer. A typical spectrum has a 1 cm ' resolution and is av¬ 
eraged over a few tens of interferograms. 

The EORMOLISM set-up was used to perform experiments 
in the monolayer and sub-monolayer regime. The experimental 
set-up consists of an ultra high vacuum (UHV) chamber (base 
pressure ~ 1 x 10 “*° mbar), containing a previously oxidised 
graphite HOPG sample (7 - 400 K, controlled by a closed-cycle 
He cryostat). Molecules are dosed onto the surface via two triply 
differentially pumped beam lines. Desorption of molecules from 
the surface is monitored using a quadrupole mass spectrometer 
(QMS, Hiden HAL-3F), positioned directly in front of the sur¬ 
face. 

Isocyanic acid was prepared in the gas phase from cyanuric 
acid (HNC 0)3 via thermal decomposition of the commercially 
avail able trimer (Aldrich, 98 %) at 650° C under primary vac¬ 
uum dRaunier et al.ll200'3^ . Small quantities of CO 2 and traces 
of CO are always present in the HNCO as a residual of the syn¬ 
thesis method. 

In the multilayer experiments presented here, HNCO was 
dosed onto the surface held at 17 K via an injection line. The 
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Fig. 1 . The infrared absorption spectra of pure HNCO bom¬ 
barded with H. The spectra are as follows: a) pure HNCO de¬ 
posited at 17 K, b) HNCO bombarded with H atoms for 140 
minutes, and c) the difference spectrum of traces a and b. Trace 
d is a reference spectrum discussed in full in the text. 


HNCO ice was bombarded with H atoms (at approximately 
300 K, with a flux of ~ lO*'* cm“^s“') produced in a molecu¬ 
lar hydrogen plasma generated by a 2.45 GHz microwave dis- 
c harge. The plasma s ource and its calibration are fully detailed 
in iTheule et al.l (1201 Ih . The hydrogenation was monitored at reg¬ 
ular time intervals using IR spectroscopy. In monolayer and sub¬ 
monolayer experiments, HNCO was dosed onto the surface held 
at 90 K via a molecular beam. This ensured that the deposited 
HNCO was not mixed with the byproduct CO 2 . The HNCO was 
cooled to ~ 10 K, then bombarded with H or D atoms (at ap¬ 
proximately 300 K, with a flux of 1 + 0.3 x 10*^ cm“^s”') pro¬ 
duced in a hydrogen or deuterium plasma in the second molecu¬ 
lar beam. After bombardment, the ices were probed using mass 
spectrometry during temperature programmed desorption (TPD) 
experiments. 


3. Results and discussion 

3.1. The multilayer, bulk HNCO regime. 

3.1.1. H bombardment of HNCO 

The spectrum of pure multilayer HNCO deposited at 17 K is pre¬ 
sented in Fig. [T] curve a. The molecular species HNCO is iden¬ 
tified via its characteristic absorption bands, as listed in Table[T] 
The most intense bands are the N-H stretching mode absorp¬ 
tions at 3554, 3362, and 3231 cm”' and the N=C=0 asymmetric 
stretching mode absorption at 2240 cm”'. A minor CO 2 contam¬ 
ination is identified via the peak at 2344 cm”'. 

The pure solid HNCO at 17 K was bombarded with H 
atoms for a total of 140 minutes (corresponding to a dose of 
~ 8 X 10'^ cm”^). The temperature of the deposition and H bom¬ 
bardment was chosen to be low enough to allow the H atoms to 
have a relatively long residency time on the surface, while being 
high enough to allow the H atoms to have a high mobility on the 
surface, to penetrate the bulk ice as deeply as possible, and to 
aid in overcoming any potential activation barrier to the hydro¬ 
genation of HNCO. The spectrum of the H-bombarded HNCO is 
presented in Figure [1] curve b, while the difference spectrum of 
pure and H-bombarded HNCO is presented in Figure [T] curve c. 
The positions of the absorption bands of the observed products 
are given in Table [1] 

There are a number of crucial points to note with regard 
to the difference spectrum of H-bombarded HNCO (Fig. [1] 
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Table 1. Fundamental infrared band positions (cm '), for the species identified in the ice during the bombardment of 
where v represents a stretching vibration, and 6 a bending vibration. 


Species 

vN-H 

vC-H 

yN=C=0„, 

0 

III 

u 

0 

II 

u 

> 

^N-Hj 

HNCO 

3554/3362/3231 

- 

2240 

- 

- 

- 

NH 2 CHO 

3313/3169" 

2886" 

- 

- 

1685* 

1385" 

OCN-/* 

- 

- 

2163“ 

_ 

- 

- 

CO 

- 

- 

- 

2136* 

- 

- 

CO 2 

- 

- 

- 

- 

2344" 

- 


“ Positively identified in the H-bombarded HNCO ice. * Tentatively identified in the H-bombarded HNCO ice. Not identified in the H-bombarded 
HNCO ice. 


curve c). First, in the wavelength region 2200 - 2300 cm“' 
the peak absorption at 2240 cm * has diminished, indica¬ 
tive of a decrease of approximately 1 to 2 monolayers (of 
580 total monolayers, c alculated assuming a band str ength of 
7.8 X 10 cmmolec ' (Ivan Broekhuizen et al. II^OiMlH . Bands 
appearing at ~ 3300 cm ' and 1680 cm ' are attributed to H 2 O, 
and the band at 2342 cm“' is attributed to CO 2 ; both of these 
species are contaminants introduced via the H plasma, confirmed 
by blank experiments where the bare gold surface was bom¬ 
barded. Figure [1] curve d is a pure H 2 O spectrum at 15 K fit¬ 
ted to the 1680 cm“^ band of the difference spectrum. Thus fit¬ 
ted, H 2 O accounts for most of the main feature in the differ¬ 
ence spectrum (at ~ 3300 cm“') as well as the full 1680 cm“’ 
band. Remaining differences between the spectra of pure and IT- 
bombarded HNCO are relatively minor. 

After inspection of the difference spectrum it appears that 
the expected hydrogenation product, methanamide (formamide, 
NH 2 CHO), has not been produced in detectable quantities. 
Arguably the simplest potential product, formamide would be 
formed by the direct hydrogenation of HNCO, as in Eq. [1] The 
left panel of Figure|2]shows two alternative fits to the 1680 cm ' 
feature (curve a); curve b is the same scaled H 2 O spectrum as 
in Figure [1] while curve c is a spectrum of pure NH 2 CHO. The 
H 2 O ice feature is a qualitatively better fit to the band formed 
during the experiment, fully accounting for its width and red 
wing, while also accounting for the 3300 cm ' feature as dis¬ 
cussed above (see Figure [1] curves c and d). However, traces of 
NH 2 CHO may have been formed during H bombardment and its 
absorption features are masked by the presence of H 2 O. No other 
peaks attributable to NH 2 CHO were observed in the difference 
spectrum. 

One further absorption feature in the difference spectrum can 
be analysed: a band appearing at 2163 cm ', as illustrated in the 
right panel of Figure for H bombardment times of 20 (curve 
d), 80 (curve e), and 140 (curve f) minutes. The slight shoul¬ 
der, centred at ~ 2136 cm ', can be attributed to CO, which is 
a trace contaminent identified in the H plasma along with H 2 O 
and CO 2 . The OCN“ ion is known to have an absorption band 
in the region of 2163 cm~ ^ (Ivan Broekh uizen et aT1l2004 e.g.) 
and could form thermally (iTheule et al 2011 ) by reaction with 
the deposited H 2 O: 


HNCO + H 2 O OCN^ + 


( 2 ) 


The counterion HaO^ was not observed in the H-bombarded 
HNCO ice, but it should be note d that the band strength is 
very weak dFalk & Giguerelll957l). and much lower than th at 


of OCN (1.3 X 10 cmion ^ ( ^n Broekhuizen et al. II2004 Ih . 


so the absorption features would be vanishingly weak in our 
ice. If OCN“ is the species responsible for the band centred at 



Fig. 2. The observed products of H -H HNCO in the multilayer 
regime. In the left panel, the spectra are as follows: a) the new 
band at 1680 cm ' (magnification of the difference spectrum in 
Fig. [I};), b) a spectrum of pure H 2 O, scaled to the band at 1680 
cm“^ (magnification of Figure [T]l), and c) a spectrum of pure 
NH 2 CHO, scaled to the band at 1680 cm"^ In the right panel, 
the spectra are the difference spectra of H-bombarded HNCO 
and pure HNCO are shown for bombardment times of d) 20 
minutes, e) 80 minutes, and f) 140 minutes (magnification of 
the difference spectrum in Fig.[T];). 


2163 cm ', it is present in very small quantities: approximately 
0.1 monolayers account for the peak produced after 140 min¬ 
utes of H bombardment. However, given the quantity of H 2 O 
deposited over the course of the H bombardment, it is not un¬ 
reasonable that HNCO could react as in E quation HI as the rat e 
of reaction was found to be 26 kJmol * dTheule et al.ll201 ih . 
which could be delivered to the ice surface by room temperature 
H atoms. 

Another potential, but rather unlikely, explanation for the ab¬ 
sorption feature is that the radical species NCO is formed by 
decomposition of HNCO, probably by H abstraction: 

HNCO + H ^ NCO + H 2 . (3) 

The spontaneous decomposition of HNCO adsorbed on 
metal surfaces has b een extensively reported in the surface sci¬ 
ence literature (e.g. iKiss & SolvrnosilH983t ICelio et ^Il997t 
Ijones & Trenarvll2008l) . The NCO radical formed by this de¬ 
composition is relatively stable, being observed on the surface 
at tem peratures up to 600 K under UHV conditions (ICelio et al.l 
I 1997 I) . The position of the absorption band associated with the 
asymmetric stretching mode absorption of NCO is highly de¬ 
pendent on the environment, with a range of around 2145 - 
2305 cm * (e.g. lNemeth et al.ir2007l) . 

One final comment to make about the multilayer reactivity 
of HNCO is that the VUV irradiation of pure HNCO has pre- 
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Fig. 3. Main panel: The infrared absorption spectra of pure 
NH 2 CHO bombarded with H atoms, a) pure NH 2 CHO deposited 
at 15 K, b) NH 2 CHO after 190 minutes of H bombardment, c) 
the difference spectrum of a and b. Inset: Infrared spectra of 
NH 2 CHO deposited at 15 K then heated at 2 Kmin"\ taken at 
15 K, 115K, and 125 K. 


vious ly been shown to pr oduce formaldehyde, formamide, and 
urea (iRaunier et al.ll200^ . In this work, urea was not observed 
in the IR spectra after H bombardment of HNCO. As the N=C 
bond must be broken to form urea, this route is unlikely. The 
only clear product of H bombardment of HNCO is the species re¬ 
sponsible for the absorption band at 2163 cm"'; probably OCN". 
The other absorption features remain ambiguous. 

3.1.2. H bombardment of NH 2 CHO 

As the results for the hydrogenation of HNCO were inconclu¬ 
sive, the reactivity of the expected product, formamide, was ex¬ 
amined. The hydrogenation of formamide could lead to the pro¬ 
duction of aminomethanol according to 

NH 2 CHO + 2H ^ NH 2 CH 2 OH. (4) 

The molecule NH 2 CH 2 OH is the most saturated form of the 
NCO moiety, so is the logical end-point for the hydrogenation 
chemistry of HNCO. 

Gas phase NH 2 CHO was dosed onto the surface held at 15 K, 
then bombarded with H atoms for 190 minutes using the RING 
experimental set-up. The spectrum of pure NH 2 CHO at 15 K is 
presented in Figure [3j curve a, with the H-bombarded sample 
spectrum as curve b. The difference spectrum of pure and H- 
bombarded NH 2 CHO is presented in Figure[3j curve c. 

Pure formamide exhibits peaks at 3313, 3169, 2886, 1683, 
1385, and 1327 cm"', as seen in Figure [3] curve a. The CO 2 
molecule was present as a pollutant in the deposited NH 2 CHO, 
as shown by the absorption feature in curve a (at 2341 cm"'), and 
in the H plasma, resulting in absorptions in curves b and c. There 
is no clear decrease observed in the peak of the vC=0 absorption 
band at 1681 cm"', and thus we conclude that NH 2 CHO has not 
reacted. 

Similarly to the bombarded HNCO, the difference spectrum 
of bombarded NH 2 CHO is dominated by absorption in the range 
3600 - 3000 cm"'. As discussed above, this might have multi¬ 
ple sources, but a large contribution is due to H 2 O. The four 
small peaks on the red wing between 2961 and 2859 cm"' are 
again attributed to contamination of the H plasma due to pri¬ 
mary vacuum. Such contamination is only observed because of 
the very long irradiation times. The peak observed in the differ¬ 
ence spectrum at 1676 cm"' is interpreted as further evidence 


of the deposition of H 2 O, with a potential contribution from the 
crystallisation of the NH 2 CHO sample during H bombardment. 
We confirmed that, after heating a sample of pure NH 2 CHO at 
2 Kmin"', the crystallisation of NH 2 CHO results in the shifting 
of the vC=0 absorption band by 4 - 5 cm"' towards the blue, 
as illustrated in the inset of Figure [2 After H bombardment of 
NH 2 CHO, no formation of OCN" was observed. We conclude 
that formamide is unreactive to H under our experimental condi¬ 
tions, and aminomethanol is not formed via this route. 

3.2. The monolayer HNCO regime. 

Because of the uncertainty in interpretation of our multilayer 
HNCO data, related to the low level contamination introduced 
by the H plasma of RING, we decided to perform experiments 
in the monolayer regime using FORMOLISM. Both experimen¬ 
tal set-ups are described in § |2 but the major advantage of 
FORMOLISM for this study is a base pressure of ~ 10"® mbar 
in the last stage of the molecular jets, which reduces the level of 
contamination of the H introduced into the chamber. Working in 
the monolayer regime allows us to bombard all deposited HNCO 
with H, unlike the multilayer regime where hydrogen penetrates 
the first ~ 1-3 monolayers, and the signal of the bulk HNCO ice 
overwhelms that of product molecules in both the Fourier trans¬ 
form infrared (FTIR) and mass spectra. However, the monolayer 
regime also limits the signal-to-noise ratio of potential products 
observed using mass spectrometry and, in these experiments, we 
were unable to observe a clear signal for our HNCO reactant 
using FTIR spectroscopy. 

3.2.1. Desorption energy of HNCO 

Using FORMOLISM, we performed experiments to determine 
the deposition conditions necessary to deposit 1 ML of HNCO 
on the graphite surface; HNCO was deposited for a series of 
fixed times, and TPDs were performed after each deposition. 
The TPD spectra are shown in Figure 0] The characteristic 
monolayer deposition time was determined by visual inspection 
of the leading edges of these spectra. All HNCO depositions 
were subsequently quantified by comparison with the identified 
monolayer deposition (Figure 01 curve b). The low flux in the 
molecular beam allows us to deposit very reproducible quanti¬ 
ties of molecular species. It takes almost 13 minutes to deposit 
1 ML of HNCO, so the deposition uncertainty on coverages of 
0.5 ML and 1 ML (used in this work) is vanishingly small and 
the depositions are highly reproducible. 

The multilayer desorption energy of HNCO was calculated 
via two different methods. The rate of desorption by un it surface, 
r, can be expressed by the Pol anyi-Wigner equation (iRedheadI 
1 19621 ICarteiin~962l : lKinglll9'7^ . where the desorption rate con¬ 
stant kde.s is described in terms of an Arrhenius law, 

dN A F iRT „ 

dT p 

where A is the pre-exponential factor, jS — ^ — 12 Kmin"' is 
the heating rate, Fjes is the energy of desorption of a molecule 
from the surface (Jmol"'), R is the gas constant (JK"' mol"'), 
T is the temperature of the surface (K), N is the number 
of adsorbed molecules on the surface (molecules cm"^), and 
n is the order of the reaction. The units of A depend on 
n: molecules'"" cm"^''^^" s"'. When analysing the desorption of 
1.5 ML of HNCO (Figure 0] curve a), zeroth order desorp¬ 
tion kinetics are assumed, as is standard practise for the mul- 
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Fig. 4. Characterisation of the HNCO monolayer deposition. The 
TPD spectra (m/z 43) for HNCO depositions of a) 1.5 ML, b) 1 
ML, c) 0.75 ML, and d) 0.5 ML. The original data are plotted in 
grey, with a smoothed version overplotted in black. The Polanyi- 
Wigner fit to the leading edge of the 1.5 ML deposition (see 
iEZB is plotted as a dotted line. 
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Fig. 5. H or D bombardment of HNCO. Left panel: TPDs of a) 
0.5 ML HNCO (m/z 43, HNCO), b) 0.5 ML HNCO + 10 minutes 
H (m/z 43). Right panel: TPDs of a) 0.5 ML HNCO (m/z 43, 
HNCO), c) 0.5 ML HNCO + 10 minutes D (m/z 43), d) as for 
spectrum c, but for m/z 44 (DNCO). The original data are plotted 
in grey, with a smoothed version overplotted in black. 


tilayer desorption of bulk material, and therefore A has units of 
molecules cm“^ s“ *. 

By fixing the pre-exponential factor. A, at a value of 
10^^ molecules cm“^ s"* (assuming that the lattice vibrational 
frequency of the solid is 10^^ s * and the number of molecules 
in a monolayer is approximately 10*^ cm“^), the desorption 
energy is calculated as the only free parameter in the fit 
of Equation |5l to the leading edge of the experimental data 
(ICollings et alJl200^ . This fit is plotted as a dotted line in 
Figure |4l Using this method, the calculated desorption energy 
was £ads,HNCO = 32.9 + 1.7 kJ moL' (3957 + 2 04 K). _ 

A second fitting method, proposed by iHasegawa et al.l 
(l1992h . assumes that the pre-exponential factor in Equation (|5]l 
is a function of Fads approximated by 

. /^A^MLEads 

A - Nml-v - Nml -W —, (6) 

where M is the mass of the adsorbate molecule, and 
Nml ~ 10*^ cm“^. The advantage of this method is that the fit 
requires only one variable, Fdes. rather th an assuming or fitting 
the pre-exponenti al factor and fitting Fdes (lAcharvva et alJ2007t 
iNoble et al.l2012h . Using the second method, the derived desorp¬ 
tion energy was Fads.HNCO = 31.0 + l.hkJmoU' (3729+ 192K). 

3.2.2. H and D bombardment of HNCO 

A series of experiments were performed on the H or D bom¬ 
bardment of monolayer and sub-monolayer quantities of solid 
HNCO at 10 - 15 K. In this temperature range, H atoms do 
not form a bulk solid because of the rapid self recombina¬ 
tion on the surface. A steady state, with an H? coverage of 
~ 10 % (lAmiaud et al.ll200'^ iKristensen et al.ll2011ll due to des- 
orptio n, is reached within the first ten seconds of H bombard¬ 
ment (ICongiu et al.ll2009l) . 

Our initial experiments into the monolayer H or D bombard¬ 
ment of HNCO were carried out on short timescales. The re¬ 
sults of 10 minutes (i.e. 1.7 ML) of H or D bombardment of 
0.5 ML HNCO are shown in Figure |5] The left panel includes 
two curves: the TPD spectrum of m/z 43 (HNCO) for 0.5 ML 
pure HNCO (curve a) and the corresponding spectrum after 10 
minutes of H bombardment (curve b). The H-bombarded HNCO 
sample contains only 0.42 ML HNCO, which represents a 16 % 


decrease compared to the pure HNCO sample. No m/z 45 (for- 
mamide) was observed to desorb during the TPD experiment, 
suggesting that no reaction occurred during H bombardment. 
However, the difference of 0.08 ML HNCO could be, at least 
partly, due to the low signal-to-noise ratio of the spectra (which 
is limited by the simultaneous measurement of multiple m/z). 

The results of D bombardment, shown in the right panel of 
Figure |5] are more revealing. After D bombardment we observe 
the species DNCO in our ice sample. In Figure|5]the TPD spec¬ 
trum of m/z 43 (HNCO) for 0.5 ML pure HNCO is again plotted 
as curve a. The remaining curves, corresponding to 10 minutes 
of D bombardment (1.7 ML), are m/z 43 (HNCO, curve c) and 
m/z 44 (DNCO, curve d), respectively. In this experiment, only 
0.33 ML of HNCO remain after 10 minutes of D bombardment, 
representing a decrease of 34 %. Assuming that the mass spec¬ 
trometer is equally sensitive to HNCO and DNCO, we calculate 
the quantity of DNCO on the surface after D bombardment to 
be ~ 0.18 ML. Thus, the total quantity of HNCO and DNCO 
on the surface represents 0.51 ML, i.e. equivalent to the quantity 
of HNCO originally deposited. We will discuss the implications 
of the deuteration of HNCO in § IT2A1 No m/z 47 (NHDCDO, 
formamide) desorbed during the TPD, therefore we observe a 
conversion efficiency from HNCO to DNCO of approximately 
100 %. 

3.2.3. The search for formamide 

In order to favour the formation of formamide, we performed 
an experiment in which we bombarded 1 ML of HNCO with 
D for 150 minutes (25 ML), comparable to the multilayer bom¬ 
bardment times (see § 13.1.11 1. The results of this experiment are 
presented in Figure |6] A reference TPD of 1 ML pure HNCO 
(m/z 43) is shown in curve a. The HNCO (m/z 43) and DNCO 
(m/z 44) on the surface after 25 ML (150 minutes) of D bom¬ 
bardment are shown in curves b and c, accounting for 0.43 ML 
and 0.22 ML, respectively. Thus, 0.35 ML (35 %) of the de¬ 
posited HNCO is unaccounted for after D bombardment. 

In order to determine the origin of the HNCO loss, we 
searched for the presence of formamide in the D-bombarded 
HNCO sam ple. The TP D method is typically sensitive to 
~ 0.01 ML dNoble et al.ll201ll) . In Figure |7] curve a, we show 
the TPD spectrum of 1 ML pure HNCO (m/z 43) as a reference. 
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Fig. 6. Deuteration of HNCO. The spectra are a) reference TPD 
of 1 ML HNCO (m/z 43), b) TPD after 25 ML (150 minutes) 
of D bombardment of 1 ML HNCO (m/z 43, HNCO), c) as for 
spectrum b, but for m/z 44 (DNCO). The original data are plotted 
in grey, with a smoothed version overplotted in black. 



Temperature / K 


Fig. 7. Presence of m/z 47 in monolayer experiments. The spec¬ 
tra are a) m/z 43 desorbing during the TPD of 1 ML HNCO, b) 
the corresponding desorption of m/z 47, c) m/z 47 after 120 min¬ 
utes of D bombardment of the graphite surface, d) m/z 47 in an 
HNCO sample bombarded with 25 ML (150 minutes) of D. 


All other curves are m/z 47 (NHDCDO). Curve b is from the 
same experiment as curve a, and indicates that there is no m/z 47 
in the deposited HNCO. Curve c is a TPD performed after 120 
minutes of D bombardment of the graphite surface, and shows 
that little or no m/z 47 is present in the atomic beam of D. Curve 
d is the TPD of the HNCO sample bombarded with D for 150 
minutes (see Figure 1^. There is no signal from m/z 47 after D 
bombardment of HNCO; we thus conclude that formamide is not 
formed at detectable levels (1 % of the deposited HNCO) during 
this experiment. 

3.2.4. Proposed reaction mechanism 

We conclude, based on the evidence from both multilayer and 
monolayer H and D bombardment of HNCO, that the expected 
hydrogenation of HNCO to NH 2 CHO (Equation [T]l does not 
readily occur and thus has a very high barrier. However, it is 
clear that D bombardment of HNCO produces DNCO and, ad¬ 
ditionally, that for longer bombardment times there is some loss 
of HNCO that is not accounted for by this DNCO formation. 
We provide a list of typical experiments, including the percent¬ 
age loss of HNCO, in Table |2] By comparing the loss of initial 
HNCO to the doses of H and D, we see that there is a general 
trend in the data, whereby one molecule of HNCO is lost from 


the sample for approximately 1 - 1.5 % of the incident H and D 
atoms. Although there was fluctuation in the absolute quantity of 
HNCO lost from the sample in our experiments, a loss was con¬ 
sistently seen for all H and D exposures greater than ~ 2 ML. 
Thes e losses are likely d ue to chemical desorption from the sur¬ 
face dPulieu et al.ll 20 i:^ . 

Chemical desorption is the stimulated desorption of reac¬ 
tant or product molecules due to an exothermic surface reaction. 
Recent experimental studies using FORMOLISM have high¬ 
lighted the significance of the chemical des orption mechanism 
in the monolayer regime at low temperat ures dPulieu et alj2013] ; 
iMinissale et alJl201^lNoble et al.ll201 Ih . It has been shown that 
up to 90 % of molecules (on a graphite surface) formed by 
radical-molecule or radical-radical surface reactions sublimate 
from the surface after formation because of the inability of 
the surface to quench the energy released during the reaction. 
On ice surfaces, more astrophysically relevant for the study of 
molecules such as HNCO (which will not be present on bare 
grains), this number is typically much lower. 

To account for both the formation of DNCO (after D bom¬ 
bardment) and the loss of HNCO from the surface (after both H 
and D bombardment), the chemistry is required to proceed via a 
cyclic pathway, such as the two proposed here; 

HNCO H 2 NCO HNCO + H 2 , (7a) 

HNCO OCN + H2^ HNCO. (7b) 

Cycle |2h corresponds to a classical hydrogenation pathway, 
i.e. H addition, but at the second step the branching ratio be¬ 
tween the products NH 2 CHO (Equation [T]i and HNCO is fully 
dominated by the reverse reaction to the initial HNCO popu¬ 
lation. Cycle |7]r corresponds to the abstraction of a hydrogen 
atom by the incident H atom. The reaction intermediates H 2 NCO 
and/or OCN are not expected to be observed in the TPDs as 
these radicals would either react with incident H atoms at a 
much faster rate than HNCO or, ultimately, recombine during 
the heating phase. It is also possible that, if cycle |7]r occurs, 
the second step could result in the isomerisation of HNCO to 
HOCN (as previousl y observed during the heating of mixtures 
of HNCO and H 2 O dTheule et alJl201 ih l. However, under cur¬ 
rent sub-monolayer experimental conditions it is not possible 
to differentiate between HNCO and HOCN, and no evidence of 
HOCN was observed in the multilayer experiments. 

The results of our D bombardment experiments are critical to 
explaining the reaction mechanism. At a low dose of D (Table |2] 
and Figure |5]l, we do not observe any loss of reactants or prod¬ 
ucts from the surface, but deuteration appears to be very efficient. 
Around 34 % of the initial HNCO (0.17 ML) is transformed into 
DNCO. Of the total D atoms incident on the surface, approx¬ 
imately 10 % have been included in the product. This repre¬ 
sents a high efficiency, particularly when we consider that the 
reaction D -H HNCO is in competition wit h D -H D. For com¬ 
pa rison, the CO + H system described by iFuchs et^ (l2009h 
or lWatanabe & Kouchil (l2002h uses only 1 % of the incident H 
atoms in the formation of H 2 CO and CH 3 OH. Although we per¬ 
form our experiments under different conditions, we are able to 
place the reaction HNCO - 1 - H between that of H - 1 - H 2 CO (which 
is faster) and H - 1 - CO (which is slower). It is clear that the barrier 
to the reaction HNCO -H D is not very high, and is probably in 
the region of 250 - 700 K. A specific study investigating varying 
doses and surface temperatures would be required to calculate 
this value; this represents a large amount of work and is beyond 
the scope of this paper. 
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Table 2. The results of experiments carried out on monolayer and sub-monolayer quantities of HNCO. 


Initial HNCO 
(ML) 

Exposed to 

Temperature 

(K) 

HNCO reduction 
(%) 

DNCO produced 
(%) 

Loss (to CD) 
(%) 

0.5 

1.7 MLD 

10 

34 ±5 

36 ±5 

-2 ±5 

1.0 

lOMLH 

7 

11 ±5 

none 

11 ± 5 

1.0 

25 MLD 

15 

57 ±5 

22 ±5 

35 ± 5 


When we extend D bombardment to longer times, we ob¬ 
serve a loss of material from the surface in addition to the for¬ 
mation of DNCO (Table|2]and Figure| 6 ]l. For an initial deposition 
of 1 ML HNCO subjected to 25 ML of D, we recover 43 % of the 
HNCO and produce 22 % DNCO i.e. a loss of 35 % of the orig¬ 
inal material. Assuming that the proportion of D atoms reacting 
is the same as for short timescales (~ 10 %), we approximate 
that HNCO + D occurs 2.5 x lO'^cm^^ times. The observed 
loss represents 3.5 x 10'^ molec cm“^, and we can make a rough 
estimation that the efficiency of c hemical desorption is ~ 14 %. 
Unlike in some previous studies (iDulieu et al.ll2013l) . the rela¬ 
tive efficiency of the desorption mechanism must be low, as only 
0.35 ML is lost from the surface over a bombardment period of 
150 minutes. This is equivalent to a desorption rate of ~ 4 xl0“^ 
ML s ', and as such any desorbing molecules would not be ob¬ 
servable using the mass spectrometer. The result indicates that 
the reaction driving the chemical desorption mechanism is not 
highly exothermic. 

Each of the two reaction scenarios that we proposed above 
proceeds via two steps. The first step (H-abstraction or H- 
addition) limits the kinetics of the overall reaction and is also 
less likely to provoke chemical desorption since less energy will 
be released than in the second radical-radical step. Thus most 
of the desorption should occur at the second reaction step. If 
the reactivity of both the HNCO and DNCO species was equiv¬ 
alent, we would assume that, with the longer exposures, over 
90 % of the initial HNCO should be transformed, which is not 
the case. The low deuteration yield could be due to the pres¬ 
ence of a second reactant (DNCO) and/or steric limitation of 
the HNCO deposition. However, we propose that the most likely 
limiting factor is the stabilisation of the HNCO by its forma¬ 
tion of dimers or larger polymers, preventing its reaction with H. 
Experiments performed on HNCO in a H 2 O ice, or isolated in a 
low temperature matricial gas, would reduce the intermolecular 
bonding and may allow hydrogenation to occur. Because of the 
inability to observe desorbing HNCO/DNCO (owing to the low 
rate of chemical desorption) and the hypothesised intermediates 
(OCN and/or H 2 NCO) by mass spectrometry, we are unable to 
determine whether one reaction cycle dominates the reactivity of 
HNCO -H H. In order to differentiate between the abstraction and 
addition mechanisms, a full ah initio quantum mechanical cal¬ 
culation treatment including surface effects would be required. 

By comparing our monolayer experiments with those in bulk 
ice, we confirm that formamide is not produced in detectable 
quantities after H bombardment of HNCO. It is also plausible 
that the inconsistency between the loss of HNCO (~ 1 - 2 ML) 
and the only confirmed product OCN /OCN (~ 0.1 ML) in our 
bulk experiments can be explained by the proposed cyclic re¬ 
activity of HNCO combined with chemical desorption. The ef¬ 
ficient H/D cyclic substitution suggests that the N-H bond in 
HNCO is relatively weak, which strengthens the argument for 
a low level of HNCO decomposition into OCN. We show our 
proposed overall reaction scheme in Figure [ 8 ] 


+H 

HNCO,,, NH,C0,„ 

HNCO,^, 


+H 

NHjCHO,,, 


Fig. 8 . Proposed reaction scheme for HNCO with H. 


4. Astrophysical implications 

Our results contradict theoretical studies which proposed that 
the initial hydrogenation step to form either HNCHO or 
NH 2 CO is very important to HNCO chemi stry and that subse¬ 
quent hydrogenation leads to NH 2 CHO (e.g. lGarrod et al.|[2008l : 
[Tides well et al.l2010h . However, ab initio calculations determine 
an activation barrier of 1390 K to the first hydrogenation step 
(iNguven et al.lll99^ . A recent experimental study of the forma¬ 
tion of HNCO by addition of the radical NH to CO revealed no 
formation of NH 2 CHO aft er co-deposition of N, H, and CO, de¬ 
spite formation of HNCO (iFedoseev et al.ll2015f) . 

The NH 2 CHO molecule has bee n tentatively identifi ed in 
ices towards the objects W 33A (ISchutte et al.l 1 19991) and 
NGC 7538 IRS9 ba sed on comparison with laboratory spectra 
(iRaunier et al.ll200^ . but the identification of molecules in the 6 
- 8 yum region is complicated by the overlap of multiple absorp- 
tion bands, and the re sults of these two studies are not conclusive 
(iBoogert et al.l20d^ . What are the reactions that would drive the 
chemistry of NH 2 CHO in the ISM? Although gas phase forma¬ 
tion routes to complex molecules exist, it is generally accepted 
that grain surface chem istry is the dominant f ormation mech¬ 
anism for such species (iBisschop et al.ll2007ah . The NH 2 CHO 
molecule can form by recombination between the radicals NH 2 
and HCO foll owing energetic processing of ices by, for example, 
UV p hotons ( Allamandola et al. 1999inMunoz Caro & Schutta 
I 2 OQ 3 I: l^n Broekhuizen et al. 112004 ). 

It has recently been determined that NH 2 CHO is the most 
energetically stable CH 3 NO isomer that can be fo rmed, and that 
the amide bond is the mos t stable bond possible (iLattelais et al.l 
12010 I) . E ones et all (1201 ll) contend that the hydrogenation of 
HNCO is insignificant as a route to formamide formation 
because of the low barrier to the reacti on HNCO -H NH 3 
(IRaunier et al.l200^lMispelaer et al.l2012l) . While it is true that 
NH 3 does react rapidly with HNCO, the flux of hydrogen atoms 
onto a grain will be superior to the quantity of NH 3 molecules 
in the ice mantle. This is particularly true before formation of 
NH 3 in the mantle. In order to reproduce observed abundances of 
HNCO and OCN' in dense molecular clouds, it has been specu¬ 
lated that hyd rogenation of HNC O must dominate its destruction 
mechanisms (iTheule et al.11^1 ll) . However, our results suggest 
that hydrogenation is not an efficient process and thus thermal 
reactions, such as the reaction with NH 3 or H 2 O, should domi¬ 
nate HNCO grain surface chemistry. 
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It has already been demonstrated that CH 3 CHO undergoes 
hydrogenation after H bombardment, but th at it also forms CH 4 , 
H 2 CO, and CH 3 OH (iBisschop et alj|^07bl) . The hydrogenation 
product, ethanol, represents approximately 20 % of the total 
products formed. The relative strength of the C=0 bond pre¬ 
vents hydrogenation being the most prevalent reaction, and the 
C-C bond breaks to form products that are chemically simpler 
and of lower mass than the original CH 3 CHO. Our results con¬ 
firm that hydrogenation of the C =0 bond is not favourable under 
low temperature H bombardment conditions. 

Our results show that the OCN moiety is formed in its rad¬ 
ical or anionic form during H bombardment. The OCN~ ion 
has been positively identified in interstellar ices (ISoifer et alJ 
Il979t iGrim & Greenber3ll987 ) and i ts rotational spectrum has 
been measured (iLattanzi et al. 1201 Oh . so searches in the gas 
phase are possible. The neutral OCN dKawaguchi et al.lll985h 
has not yet been identified in th e ISM; it has, however, l ong 
been predicted in dense clouds (iPrasad & Huntres^ Il978h . If 
OCN/OCN“ can be formed by the reaction of H with HNCO, 
this could help explain the destruction of HNCO, while con¬ 
firming that formamide formation does not occur via this route. 
Additionally, in dense regions such as pre-stellar cores, H/D ex¬ 
change is potentially the do minant destruction mechanism for 
HNCO (iRoberts et alJl200^ . 

We have also shown that formamide does not react with hy¬ 
drogen to produce aminomethanol, the logical saturated end¬ 
point for HNCO hydrogenation. Aminomethanol has, however, 
already been shown to form under interstellar conditions by the 
purely thermal reaction between H 2 CO and NH 3 , with an acti¬ 
vation barrier of 4.5 + 0.5 kJmol ' (541 + 60 K, iBossa et ^ 
(l2009lH . Theoretical studies suggest that, in the presence of 
acids, aminomethanol can undergo spontaneous dehydration 
to fo rm methylenimine (CH 2 =NH) and water (IWalch et alJ 
l 200 lh : experiments have shown that aminomethanol yields 
hexamethylenetetramine (HMT) from t he polymerisation of 
the dehydration produc t methylenimine jBernstein et al.l fl 995t 
IVinogradoff et alJl 201 ih . Although it has not yet been observed 
in the interstellar medium, in the la boratory HMT has a sub- 
limation temperature of ~ 553 K (iBernstein et alJ 1 19951) . so 
it is thought to be present as part of an organic residue on 
comets. The formation routes to aminomethanol and hydrox- 
yacetonitrile (HOCH 2 C N) are competitive when CN“ is in¬ 
cluded in the initial ice dPanger et al.ll20l3) . These destruction 
routes and competitive reactions, combined with our conclu¬ 
sion that aminomethanol does not form by the hydrogenation 
of formamide, might result in low interstellar abundances of 
aminomethanol, although it has not yet been observed. 
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